Dynamics of matter solitons in weakly modulated optical lattices 
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It is shown that matter solitons can be effectively managed by means of smooth variations of 
parameters of optical lattices in which the condensate is loaded. The phenomenon is based on the 
effect of lattice modulations on the carrier wave transporting the soliton and that is why is well 
understood in terms of the effective mass approach, where a particular spatial configuration of the 
band structure is of primary importance. Linear, parabolic, and spatially localized modulations 
are considered as the CEise examples. It is shown that these defects can originate accelerating and 
oscillating motion of matter solitons as well as simulate soliton interaction with attractive and 
repulsive defects. 

PACS numbers: 



X 
J3 



I. INTRODUCTION 



Realization of a Bose-Einstein condensate (EEC) in an 
optical lattice 1] originated intensive experimental and 
theoretical studies Q,lli,lllllllllili[l3,[ll|ll 

IT^ IT^ of the phenomenon. One of the main features 
induced by the periodicity is the appearance of a band 
structure in the spectrum of the underlying linear sys- 
tem, i.e. the gas of noninteracting atoms. The band 
spectrum is responsible for a number of effects, including 
modulation instability's','?! ^^^d formation of gap solitons 
, Landau-Zener tunneling |^ , Bloch oscillations of 
EEC's 0,|EI3' lensing effect soliton stabilization 0, 
etc. These effects, being based on the properties of the 
linear system, are well described in terms of the concept 
of the effective mass, which takes into account wave na- 
ture of the phenomenon and that is why its inverse value 
is also referred to as the group velocity dispersion 4] . In 
practice, laser lattices are never perfect. In particular, 
they are often imposed simultaneously with other exter- 
nal potentials, like, for example, a magnetic trap which 
significantly affects the motion of a soliton • Moreover 
one can create spatially localized lattices, propagation of 
matter waves through which displays a number of inter- 
esting properties such as resonant transmission and 



soliton generation through the modulational instability. 
It is also highly relevant to mention that very recently 
direct experimental observation of a gap matter soliton 
in ^^Rb EEC, containing about 250 atoms, has been re- 
ported [il]. 

In this paper we present analysis of the EEC dynamics 
in an optical lattice with smoothly modulated parame- 
ters. In this case the band gap structure is still pre- 



served but is deformed by the modulation. In particular, 
we show that this is a way of managing the dynamics 
of matter waves - making them accelerating, oscillating, 
etc. 

The model and the physics of the phenomenon are de- 
scribed in Sec|nJ Acceleration of a matter soliton in a 
linearly modulated lattice is considered in Sec. IIIII Soli- 
ton oscillations in a lattice with a parabolic modulation 
of the depth are considered in Sec. II VI Sec. is devoted 
to soliton interaction with a spatially localized defects. 
The results of the paper are summarized in the Conclu- 
sion. 



II. PHYSICS OF THE PHENOMENON 



Let us consider a trap potential which can be written 



'+Loly^+i^lz^) + V^{x). 
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down in the form Vtrap ~ ^{'^\\^' 
The first term in the right hand side of this expression 
describes a magnetic trap with t^n and being the longi- 
tudinal and transverse linear oscillator frequencies. The 
condensate is chosen to have a cigar shape with the as- 
pect ratio satisfying the relation wh/wq ^ Oo/S^ ^ 1 ('^o 
and ^ being the transverse linear oscillator length and 
the heahng length, respectively). We are interested in 
excitations of a EEC having characteristic scales of order 
of the healing length and a relatively small amplitude, 
what allows us to neglect the term yw^x^ in the expres- 
sion for the trap potential Vtrap- The term Ve(a;), where 
e is a deformation parameter controlling the potential 
shape, describes smoothly modulated optical lattice. It 
will be assumed to have the form V^{x) = /(e'^/^a;)Vo(a;), 
where Vq{x) = Vq{x -\- L) is an unperturbed lattice hav- 
ing a period L of order of the linear oscillator length cq, 
L ~ flo, and f{e^^^x) is a smooth modulation normalized 
as follows: /(O) = 1. Smoothness in the present context 
means slow variation of f{e^^^x) on the scale of the heal- 
ing length, or in other words e ~ oq/C- Then, one can 
apply the multiple-scale expansion [j| in order to reduce 
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the Gross-Pitaevskii equation governing the dynamics of 
the BEG to an effectively one-dimensional (ID) nonlin- 
ear Schrodinger (NLS) equation with the slowly varying 
effective mass, rua = ma{e^^'^X) = [dl£a{k;e^^^X)]^^, 
and the group velocity of the carrier wave depending in 
the coordinate, Va = dkSa{k] e^^'^X): 



(1) 



Hereafter a and k stand for a number of the zone and for 
the dimensionless wave vector in the first Brillouin zone 
(BZ), respectively, x =signas, Os being the s-wave scat- 
tering length. In ^ we have passed to dimensionless slow 
variables: X — tx/ao = ei and T = e^uj^t = e^t. Tak- 
ing into account smoothness of the potential, the spec- 
trum, £a{k; e^^^x) in some point of the space, say x = xq 
(xq = Xo/ao), can be computed from the linear eigen- 
value problem 



(P^ak 

dx^ 



f„(fc; - /(e3/2io)V^o(i)l ip^k = 0. (2) 



In order to get qualitative picture of effects which can 
be observed in the matter wave dynamics we carry out 
numerical simulations using a ID model (justified for for 
the low-density EEC's) 



-V'i£ + /(e'/'i)V^,(i)V + x|^|V. (3) 



It is to be emphasized here that -0 in Eq. Q describes 
the complete macroscopic wave function whereas ^ in 
Eq. J^l describes behavior of its envelope. 

In the case of a homogeneous lattice, i.e. Eq. Q with 
f{x) = 1, there exists a solitary wave solution (bright 
matter soliton) which can be found by means of the sta- 
tionary ansatz ilj{x,t) = M(i)e~*^*, where u{x) is a real 
valued function satisfying the equation as follows 



"Sx + [£ - Vo{x)] u - xu^ = 0. 



(4) 



Let us consider now the boundary of the BZ, i.e. 
|fco| = TTtto/L. As it is known (see e.g. stationary 
solitary wave solutions subject to zero boundary condi- 
tions lim|j|_,oo M(i, t) = 0, exist only if £ belongs to a 
gap of the spectrum of For the sake of dcfiniteness, 
below we are dealing with the first lowest forbidden gap, 
which lower and upper edges will be denoted as S'--^"^ and 
f^^-*, respectively (see an example in Fig. ^). Thus the 
frequency of a static matter soliton must satisfy the con- 
dition < £ < £'-^\ Then in the case x = -1 (x = 1) 
a small amplitude bright matter soliton can be excited 
with energies inside the gap in the vicinity of the upper, 
£^^^ (lower, S^^-*) band edges, correspondingly 4]. For the 
next consideration it is important to mention that enve- 
lope solitons can also exist at f > f^^-' if X = s^nd at 
£ < f^^-* if X = 1- Respective solitons are created against 
a moving carrier wave background (with |fc| |fco| and 
Va 7^ 0), i.e. they move with the velocity Uq,, and are 
described by the formula 



^s{X,T) = 



V-XmgjVexp (|m„iV^r) 
2 cosh [^rriaN {X ~ Xo{T))] 



(5) 



with either x = ^1, a = 2, and TO2 > 0, or x = 1, = 1, 
and mi < 0. Xq{T) = —VaT is a coordinate of the soliton 
center and is a normalized number of atoms 



N = 



l^l^dX. 



(6) 



Strictly speaking Eq. is a solution of the unper- 
turbed NLS equation |^ with Va and TOq constants. 
Since, however, the parameters are changing slowly in 
space, the simplified picture based on the analytic form 
(O appears to be good enough for qualitative under- 
standing various phenomena. 

In what follows the consideration is restricted to the 
case X = — 1 and discuss the dynamics of small- amplitude 
bright matter soliton near the upper edge of the first band 
gap £^'^\ 



III. MATTER WAVE ACCELERATION. 

We first consider linear modulation of the lattice am- 
plitude (see Fig. ^) 

= (l-e3/2^)cos(2i). (7) 

which lowest forbidden gap is shown on Fig. ^i. 
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FIG. 1: (a) The linearly modulated periodic potential Ve given 
by (|7|l. (b) The initial profile of the condensate near the 
upper bound of the gap for E^'^^ = 1.46. (c) The profile of 
the condensate at time t = 100. (d) The band structure the 
potential (O (solid line) and the energy of the soliton (dashed 
line). All data are computed for e = 0.02. 

Let us assume that in the vicinity of the origin, i.e. 
near a: = 0, a bright static gap soliton with i;2 = is 
created (see Fig. QJa). Due to nonzero extension of the 
soliton, it partially occupies the space where the soliton 
energy falls into allowed band (see Fig. ^) and thus cor- 
responds to running linear waves. This is a reason for a 
soliton to start to move toward the region with a narrower 
gap, i.e. in the positive direction. Since it is assumed that 
in the leading approximation the energy of the soliton is 
conserved, the group velocity of the background increases 
and, therefore, soliton velocity grows when the coordi- 
nate of the soliton center increases (Fig. [2^), what gives 
rise to a soliton acceleration (see Fig. |2l3). An example 
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FIG. 2: (a) The coordinate of the soliton center of mass, (i), 
vs time, (b) Dependences of the carrier wave group velocity 
U2 obtained directly from the periodic structure (dots) and 
velocity of the soliton center of mass obtained from the direct 
numerical simulations (solid line), (c) Evolution of the dis- 
persion of the soliton. (d) Dependence of the effective mass 
m2 on the soliton coordinate. Curves 1, 2 and 3 correspond 
to e = 0.02,0.03, and 0.04. 



of dynamics of the coordinate of the condensate center 

1/2 

of mass, (i), and its dispersion a = (^(i^) — (i)^^ , 

obtained by numerical integration of Eq. © , are shown 
in Fig. 12 (a) and (c) (angular brackets stand for spatial 
average). From this figure one can see that the soliton 
indeed undergoes acceleration in such a way that the soli- 
ton velocity follows the change of the group velocity of the 
background. Simultaneously, the dispersion of the soli- 
ton increases as it is shown in Fig. [21:. From the first sight 
this result does not resemble a model given by Eq. © 
because according to this model the soliton width should 
decrease as the effective mass increases. 

In order to explain the apparent discrepancy let us 
consider the dispersion relation of the second zone, Ei(k\ 
in the vicinity of the band edge in more detail. Since 
is an even function of the wave vector, one can 
expand 



(8) 



where f'^' = £2(^0)- We assume that (i) the center of 
soliton is displaced from the point 5 = 0, where its en- 
ergy was £27 to some point i, = X without change of the 
energy, and that ii) the gap is large enough and within 
the range < i < X in the leading order one can ne- 
glect the change of the functional dependence of Ei on 
the wave vector k. Then one has the situation illustrated 
in Fig. O When soliton moves from i = to i = X, 
the acquired energy shift, , toward the allowed zone 
can be estimated as where /3 is a small (/3 ^ 1) 
angle between the slope of the band-gap edge and the x- 
axis corresponding to the energy Ei (see the left figure in 



Fig.O. The same energy shift on the diagram £2(^)7 giv- 
ing the energy f '^^^ -I- A£, corresponds to the wave vector 
fci, which can be approximately found from the equa- 



tion (fci — fco)^ = 2mi^E: k\ = ko 



2m2l3X. Hence, 



one computes that dependence of the group velocity on 
a coordinate x is given by 



V2 



S:>0 
0, i < 



(9) 



while the dependence of the effective mass is approxi- 
mated by 



TO2(i) = m2o 



m2Qi^Am, Am 



m2odt£2{ko){l0) 



where m2o is the effective mass at the boundary of the 
zone, i.e. at i = (notice that it follows from Fig. |2l 
dtE2{ko)<Q). 

It is remarkable that the above simplified estimates 
give rather precise estimates for the relevant quantities. 
Indeed, for e = 0.02 from Fig.^ and Fig.[51i one obtains 
(3 ~ 0.0013 and m2 ~ 0.1 what gives for X = 30 an 
estimate V2{X = 30) « 0.88, in accordance with jnj while 
the numerical value is V2{X = 30) w 0.87 (see Fig. Ela). 
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FIG. 3: Schematic illustration of the calculus of the depen- 
dence of the group velocity and the effective mass on the 
coordinate, for linearly modulated potential. The notations 
are defined in the text. 

Returning to the soliton dynamics shown in Figs. ^ 
13 let us neglect small change of the effective mass, i.e. 
make a substitution — rriao in the equation for slowly 
varying amplitude and compute the following inte- 
gral identities 



MxlM^dX; 



dN 

dP — 
— = 0, P=-i I ^^xdX; 



(11) 
(12) 



d{X^ 



dT m 



— I X(«'x* - *^'x)d^ 



{X\^)xm^dX, {X^) = / X^l^ppdX 

-00 J —00 

(13) 

where N is the number of particles given by ©, P is the 
wave momentum in the frame moving with the group ve- 
locity of the carrier wave, and ^/ (X^) /N is an average 
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width of the wave packet. Then it follows from H13|l that 
the total momentum is preserved. Moreover for a soli- 
ton ansatz © it is zero, what reflects the fact that the 
condensate is moving with the group velocity of the car- 
rier wave what corroborates with the findings shown in 
Fig. Et). From Hll|l it follows that the number of parti- 
cles is not preserved. This, however does not contradict 
the conservation law of the total number of atoms in the 
Gross-Pitaevskii equation since the model given by the 
Eq. describes the main approximation of the ground 
state and does not account for the high-frequency ra- 
diation while, on the other hand, in the leading order, 
the soliton is moving in the inhomogeneous medium and 
thus represents radiating matter wave Considering 
the second band, a = 2, and approximating the solution 
by ©7 ® one can rewrite ((TT|l as 



dN 
'dT 



2m 



20 Jo 



dX 



(|*,(x,T)|2- |«',(o,r)|2) < 



what reflects the fact that the soliton looses particles. 
Finally, taking into account that 'i's, given by lO is real, 
one obtains from 



d{X^ 
dt 



^ 5 



/3 



2m2o 



/ X^'^\^s{X,T)\'^dX >Q 
Jo 



what means increase of the soliton width, observed in 
Figs. [It and 12;. 

The wave packet during its motion follows the carrier 
wave and thus rapidly acquires relatively large group ve- 
locity which is of order of unity. The latter factor does 
not allow adiabatic adjustment of the soliton parameters 
to increase of the effective mass of the carrier wave (see 
(|10|l and Fig. QJi) . Nonadiabaticity of the process mani- 
fests itself in radiative losses and broadening of the wave 
packet. This process continues as soliton moves to the 
region with small potential depth, as it is observed in 
Fig.^, and with time, as matter of fact, the wave looses 
its solitonic properties transforming itself in a dispersive 
wave packet. It worth pointing out that another simple 
way to look at the broadening of the wave packet is to 
take into account that soliton wave front has larger ve- 
locity than the velocity of its tail leading to spreading of 
the pulse. 



parabolic modulation. The upper bound of the first gap 
of this potential for e = 0.02 is depicted in Fig. 0];. 
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FIG. 4: (a) The inhomogeneous periodic potential Ve given 
by (O with e = 0.02 and Ai = 45. (b) The initial profile of 
an envelope soliton (the same as in Fig. ^) . (c) The upper 
edge of the first band gap for for the potential shown in (a). 

Now one observes oscillations of the condensate cloud 
(see Fig. EJ. The accelerating part of this motion when 
the soliton is moving toward the center of the potential 
is explained in the previous Section. The soliton having 
passed the central part of the potential is decelerating 
and at some point the velocity of the center of mass of 
the cloud becomes zero. This occurs in the vicinity of 
the turning point, Xturm where the energy of the soliton 
falls into the forbidden gap (the intersections of dashed 
and solid lines in Fig. Et). Indeed, considering a soliton 
moving from the right to the left, it follows from Eq. I®, 
that the velocity of the background, and hence the veloc- 
ity of the soliton, goes to zero as ^/x — Xturm when the 
coordinate of the soliton center approaches the turning 
point xturn- Thus one can speak about Bragg reflection 
of the soliton, from inhomogeneous periodic potential. 
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IV. OSCILLATION OF THE MATTER WAVE IN 
A LATTICE SUBJECT TO PARABOLIC 
MODULATION. 

Let us consider now a soliton dynamics (the initial 
profile is shown in Fig. 0Jd) in a lattice modulated by 
a parabolic function (see Fig. ^) 



V,{x) =\l + f-^^^{x - Ai)2j cos(2i) 



(14) 



where Ax is introduced for an initial relative shift 
between the soliton center and the minimum of the 



FIG. 5: Dynamics of the center of mass (a) and of dispersion 
(b) of the matter soliton. Parameters the same as in Fig. |1] 
The dashed line corresponds to {x} — Ax. 

There are two features of the soliton dynamics to point 
out in this case. First, in contrast to the case of the lin- 
ear modulation, at the initial stages the dispersion de- 
creases (see FigEt). This compression of the pulse is ex- 
plained by its large extension, when the leading part of 
the wave packet undergoes Bragg reflection and moves in 
the negative direction while the soliton center still moves 
in the positive direction. The second feature demon- 
strated in Fig. \^ can be identified as beatings of the 
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mean coordinate of the wave packet. In order to ex- 
plain this effect, we assume that the solution can be 
described approximately by the formula ^ where the 
coordinate of the soliton center Xq{T) depends period- 
ically on time. This dependence is determined by the 
effective linear oscillator frequency ly associated with the 
effective parabolic potential, the latter being estimated 
as = 2e5/4(« 0.014 in our simulations). This frequency 
approximately coincides with the higher frequency that 
can be seen in Fig. [SJl, and which is ^ 0.017. The ob- 
served discrepancy is due to the fact that one deals with 
a modulated periodic potential, rather than with a real 
parabolic one. On the other hand, the envelope soliton is 
characterized by the internal frequency ujint — m2N'^/8 
[see Eq. (0)]. In our numerical simulations it is a low 
frequency uji„t ~ 0.002. Thus by expanding the modu- 
lus of the soliton solution in the Fourier series with re- 
spect to time the motion is determined by superposition 
of the lowest harmonics with v ± tOint which results in 
the beating phenomenon. The above estimates imply 
'-^int/^ ~ 8.6, which corresponds to 9 fast oscillations per 
one period of the slow beatings observed in Fig. 



V. INTERACTION OF A SOLITON WITH A 
DEFECT. 

The phenomenon of the Bragg reflection of a soliton 
can be directly observed in a system where the optical 
lattice has local increase of the depth, as it is shown in 
Figs. 



1 



cos(2i). 



(15) 



Local deformations of an optical lattice will be referred 
to as lattice defects. In this before, we also start 

with a numerically obtained bright soliton, but adding 
now some initial velocity v with respect to the stationary 
background which is introduced by means of the phase 
factor V(i,0)e*"*/2^ We are interested in the scattering 
process at different initial velocities, which are assumed 
to be small enough: v <^1. 

In this case solitons are reflected from the defect (see a 
snapshot Fig|S|3). This is the effect of the Bragg reflection 
which is explained in Fig. |HJd. Indeed, due to increase of 
the lattice depth, the gap is also increased locally. That 
is why the soliton, which initially was moving with a 
relatively small group velocity u, at some point reaches 
the (curved) band edge (it is given by the intersection 
of a dashed and solid lines in Fig. ^p) which prevents 
further propagation of the carrier wave. Thus, the defect 
described by the Eq. ((T3|l can be classified as repulsive. 

Let us consider now an opposite situation, when the 
depth of the periodic potential is locally decreased (see 

Fig.EK) 



1 



cos(2i). 



(16) 
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what leads to local narrowing of the forbidden gap. 



FIG. 6: (a) The lattice defect given by Eq. lO with e = 0.05 
and Ax — 100. (b) The upper edge of the lowest gap (solid 
line) and the energy of the soliton with the initial velocity 
V — 0.1 (dashed line). Intersections of these lines represent 
the turning points (c) Profile of the reflected soliton at time 
t — 150. The initial soliton profile is the same as in Fig. 



By using the arguments based on the band gap struc- 
ture outlined below one may argue that the modulation 
described by the Eq. (|16|l acts as an attractive impurity. 
Indeed, in the region of the defect there exists a forbid- 
den band shrinking. Assuming as before, that the soliton 
energy £ is constant, one concludes that in the region of 
the defect the energy shift A£ (see Fig. toward the 
allowed zone depends on the coordinate A£ = A£{x): it 
increases as the soliton approaches the center of the de- 
fect and then decreases as coordinate increases. As far as 
larger A£(x) correspond to larger velocities of the carrier 
wave, the soliton is accelerating when it approaches the 
defect and decelerating when soliton moves outward the 
defect. 

Since the lattice modulation H16|l acts as attractive a 
number of atoms could be captured by such a defect, in 
the case when an initial kinetic energy of the conden- 
sate is small enough. This is exactly what we observe 
in numerical simulations shown in Fig. \7}p. The higher 
velocity matter waves pass through the defect without 
substantial changes (see Fig.[33). 



VI. CONCLUSION 

In conclusion, we have shown that in smoothly modu- 
lated optical lattices that can be created by using quasi- 
monochromatic laser beams, one can effectively manage 
matter solitons as they accelerate, deccelerate, oscillate 
or undergo the reflection depending on the type of the 
modulation introduced. Since the above processes are 
controlled by the periodic structure, not only dynami- 
cal, but also other properties of matter waves such as 
the energy, the effective mass and width of the soliton 
can be manipulated. Although in a particular dynamical 
process a matter wave can loose its solitonic properties, 
the effective mass approach provides a qualitative expla- 
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nation of the main features of the soliton dynamics, if 
the wave packet possesses a substantially larger exten- 
liion than the lattice period and lattice modulations are 
smooth enough. 
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FIG. 7: (a) Periodic potential with the defect given by 
Eq. (1161 and corresponding upper edge f'^' of the first band 
gap (here the dashed hne is initial soliton energy). Param- 
eters of the defect are e = 0.05 and AS — 100. Profiles of 
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initial soliton profile is the same as in Fig. |S| 
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